Signature of the long range triplet proximity effect in the density of states 
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We study the impact of the long-range spin-triplet proximity effect on the density of states (DOS) 
in planar SF1F2S Josephson junctions that consist of conventional superconductors (S) connected by 
two metallic monodomain ferromagnets (Fi and F2) with transparent interfaces. We determine the 
electronic DOS in F layers and the Josephson current for arbitrary orientation of the magnetizations 
using the solutions of Eilenberger equations in the clean limit and for a moderate disorder in ferro- 
magnets. We find that fully developed long-range proximity effect can occur in highly asymmetric 
ferromagnetic bilayer Josephson junctions with orthogonal magnetizations. The effect manifests it- 
self as an enhancement in DOS, and as a dominant second harmonic in the Josephson current-phase 
relation. Distinctive variation of DOS in ferromagnets with the angle between magnetizations is 
experimentally observable by tunneling spectroscopy. This can provide an unambiguous signature 
of the long-range spin-triplet proximity effect. 

PACS numbers: 74.45.-l-c, 74.50.-|-r 



I. INTRODUCTION 

Long-range spin-triplet superconducting correlations 
induced in heterostructures comprised of superconduc- 
tors with the usual singlet pairing and inhomogeneous 
ferromagnets have been predicted recently. ^^"^ It has been 
experimentally verified that superconducting correlations 
can propagate from SF interfaces with a penetration 
length up to 1/xm, and provide a nonvanishing Joseph- 
son supercurrent through very strong ferromagnets i^^— 

The simplest Josephson junction with inhomogeneous 
magnetization is SF1F2S structure with monodomain fer- 
romagnetic layers having noncollinear in-plane magne- 
tizations. However, the long-range Josephson effect is 
not feasible in the junctions with only two F layers,—"— 
except in highly asymmetric SF1F2S junctions at low 
temperatures, as it was shown in Refs. [l^ and [TtI . In 
that case, the long-range spin-triplet effect manifests it- 
self as a large second harmonic {I2 ^ h) in the ex- 
pansion of the Josephson current-phase relation, = 
Ii sin(0) + I2 sin(20) -I- • ■ •. The ground state in Joseph- 
son junctions with ferromagnetic barriers can be either 
or TT state. The energy of the junction is proportional to 

I{(f)')d(f)' , hence a second harmonic leads to degenerate 
ground states at (f> = and (/> = tt. Small contribution 
of the first harmonic lifts the degeneracy which results in 
coexistence of stable and metastable and tt statesj^^ 

In SFS Josephson junctions with homogeneous magne- 
tization, the projection of the total spin of a pair to the 
direction of magnetization is conserved and only spin- 
singlet and triplet correlations with zero spin projection 
occur. These correlations penetrate into the ferromag- 
net over a short distance (determined by the exchange en- 
ergy h) — hvp/h in the clean limit and ^i? ~ \JhD jh 
in the dirty limit, where diffusion coefficient D — vplj^i. 
For inhomogeneous magnetization, odd-frequency triplet 
correlations with nonzero (±1) total spin projection are 



present as well. These correlations, not suppressed by the 
exchange interaction, are long-ranged, = hvp/ksT 
{^JhD /kBT) in the clean (dirty) limit, and have a dra- 
matic impact on the Josephson effect.'^ 

Dominant infiuence of long-range triplet correlations 
on the Josephson current can be realized in SFS junctions 
with magnetically active interfaces narrow domain 
walls between S and thick F interlayers with misaligned 
magnetizations;^"— i22i"— or superconductors with spin 
orbit interaction)^ In addition to the impact on the 
Josephson current, the odd-frequency triplet pair correla- 
tions can be seen through enhanced electronic density of 
states by tunneling spectroscopy J2, It has been found that 
proximity effect in diffusive FNS or NFS structures (N is 
a normal nonmagnetic metal), with magnetically active 
interface or precessing magnetization, gives a clearcut 
signature of the odd-frequency superconducting correla- 
tionsjSiiS^i In general, measuring DOS is a very powerful 
tool to characterize the nature of superconducting corre- 
lations in NS and FS structures i^^i'^° 

In this article, we study the proximity effect and infiu- 
ence of spin-triplet superconducting correlations in clean 
and moderately diffusive SF1F2S junctions with trans- 
parent interfaces. Magnetic interlayer is composed of two 
monodomain ferromagnets with arbitrary orientation of 
in-plane magnetizations. We calculate density of states 
(DOS) in ferromagnetic layers and the Josephson current 
from the solutions of the Eilenberger equations. Previ- 
ously, F1SF2 and SF1F2S junctions with monodomain 
ferromagnetic layers having noncollinear in-plane magne- 
tizations have been studied using Bogoliubov-de Gennes 
equatio n^^i'^^" — and within the quasiclassical approxima- 
tion in diffusivei^i^"— and cleani^'"'^ limits using Usadel 
and Eilenberger equations, respectively. 

We present analytical solutions in stepwise approxi- 
mation in the clean limit, and numerical self-consistently 
obtained results both in the clean limit and for a moder- 
ate disorder in ferromagnets. The second harmonic in the 
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Josephson current-phase relation is dominant for highly 
asymmetric SF1F2S junctions composed of particularly 
thin (weak) and thick (strong) ferromagnetic layers with 
noncoUinear magnetizations. This is a manifestation of 
the long-range spin-triplet proximity effect at low tem- 
peratures, related to the phase coherent transport of two 
Cooper pairs. lS.ii2, The proximity effect is also accompa- 
nied by pronounced change of both shape and magnitude 
of electronic DOS in ferromagnetic layers near the Fermi 
surface. In this case, the self-consistency and finite (but 
large) electronic mean free path do not alter the analyti- 
cal results qualitatively. On the other hand, in symmetric 
SFFS junctions, with equal and strong magnetic influ- 
ence of F layers, the long-range proximity effect is weak 
(cannot be dominant-'-^); we show that DOS is therefore 
equal to its normal metal value even for moderately dif- 
fusive ferromagnets. 

The article is organized as follows. In Sec. II, we 
present the model and the equations that we use to cal- 
culate the electronic DOS and the Josephson current. In 
Sec. Ill, we provide results for different ferromagnetic 
layers thickness and orientation of magnetizations. The 
conclusion is given in Sec. IV. 



II. MODEL AND SOLUTIONS 

We consider a simple model of an SF1F2S hetero- 
j unction consisting of two conventional (s wave, spin- 
singlet pairing) superconductors (S) and two uniform 
monodomain ferromagnetic layers (Fi and F2) of thick- 
ness di and d2, with angle a = q;2 — ai between their 
in-plane magnetizations (see Fig. [T]). Interfaces between 
layers are fully transparent and magnetically inactive. 

We describe superconductivity in the framework of the 
Eilenberger quasiclassical theory?^>^ Ferromagnetism is 
modeled by the Stoner model, using an exchange energy 
shift 2h between the spin subbands. Disorder is charac- 
terized by the electron mean free path / = vpT, where 
T is the average time between scattering on impurities, 
and Wi? is the Fermi velocity assumed to be the same in 
S and F metals. 

Both the clean and moderately diffusive ferromagnetic 
layers are considered. In the clean limit, the mean free 
path I is larger than the two characteristic lengths: the 
ferromagnetic exchange length = hvp/h, and the su- 
perconducting coherence length ^5 = Hvf/t^^O: where 
Ao is the bulk superconducting pair potential. For mod- 
erate disorder, < I < S,s- 

In this model, the Eilenberger Green functions 

gaa'{x,0,Un), gla' , 9 , LOn) , faa'{x,6,UJn), and 

flcr'{x,0,uJn) depend on the Cooper pair center-of- 
mass coordinate x along the junction axis, angle 9 of 
the quasiclassical trajectories with respect to the x axis, 
and on the Matsubara frequencies ujn = TrkBT(2n + 1), 

n = 0, ±1, Spin indices are a =t,i- The Eilenberger 

equation in particle-hole spin space can be written in 




FIG. 1. Schematics of an SF1F2S heteroj unction. The mag- 
netization vectors lie in the y-z plane at angles qi and 02 
with respect to the z-axis. 
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with normalization condition g^ — 1. We indicate by ■•■ 
and • " • 2x2 and 4x4 matrices, respectively. The brack- 
ets (...) = (1/2) {...)sin{9)d9 denote angular averaging 
over the Fermi surface, [ , ] denotes a commutator, and 
Vx = "vp cos 9 is the projection of the Fermi velocity vec- 
tor on the X axis. 

The matrix of quasiclassical Green functions isi^^ 



(2) 



and the matrix V is given by 

y = i (g) Re h{x) a + 2x3 (g) Im h(a;) • S , (3) 

where the components crx,(Ty,&z of the vector S, and 
are the Pauli matrices in the spin and the 
particle-hole space, respectively. The in-plane {y-z) mag- 
netizations of the neighboring F layers are not collinear 
in general, and form angles ai and with respect 
to the z-axis in the left (Fi) and the right (F2) ferro- 
magnets. The exchange field in ferromagnetic layers is 
h(a;) = /ii(0, sinai, coscki) and /i2(0, sin 012, cos 012). 

We assume that the superconductors are identical, 
with 



<T2A 
(72 A* 



(4) 



for X < —di and x > d2- The self-consistency condition 
for the pair potential A = A(x) is given by 



A^~i\nNoiO)kBTj2{fn)^ 



(5) 



where A is the coupling constant, No{0) = rnkp/T^'^h^ is 
the density of states (per unit volume) of the free electron 
gas at the Fermi level Ep = h^kp/2m, and kp = mvp/h 
is the Fermi wave number. In F layers A = 0. 
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The density of states normalized by its normal-state 
value is given by analytical continuation 

N{s)/No{0) = 2 E(I^« 



We consider only transparent interfaces and use conti- 
nuity of the Green functions as the boundary conditions. 
Analytical solutions of Eq. ([T]) are derived (see Appendix) 
in the stepwise approximation for the pair potential 



where a small imaginary part of the energy is introduced 
to regularize singularities. In numerical calculations we 
choose (5/Ao(0) = IQ-^. 

The supercurrent is obtained from the normal Green 
function through the following expression 



= i7re7Vo(0)^fcBTEEKlm.g,,), (7) 



where is the macroscopic phase difference across the 
junction, and S is the area of the junction. In our ex- 
amples, the current is normalized to nAo/eRN where 
Rn = 2n^h/Se^kj,. 



i<p/2 



e(x - da) 



(8) 



where Aq is the bulk pair potential, and 0{x) is the 
Heaviside step function. The temperature dependence 
of the bulk pair potential Aq is given by Ao(T) = 

Ao(0) tanh (l.7A^Tc/T ~ 1^ M 

We have chosen a2 — 0, that is a = ai. For three 
characteristic values of the angle between magnetizations 
a = 0, 7r/2, and tt, the normal Green functions in F2 layer 
are x-independent and can be written in a compact form 



3 X UJn Al 
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Aq a;„[cos {(j) — i^n) — cos (i?i) cos (z^a)] ~ sign(cos 9)[± cos (di) sin (iJa) 



fin (f^^ -I- uj^) COS {4> - zS„) -I- (fi^ - oj^) COS (i^i) COS {^2) + 2iujnfln sign(cos 9) sin {(j) - «S„) ' 

(10) 



UJn Ag 



Wn -I- iQn sign(cos 9) cot (^±^2-^ + i- 



(11) 



respectively. Here, ± corresponds to a =t,i, fi'^ = 
(jj^ + Aq, En = 2{di + d2)uJn/fivF cos 9 and dj = 
2djhj /hvp cos 9, j = 1, 2. 

The well-known normal Green functions for a = and 
TT are the same in Fi and F2 layers, Eqs. ^ and ([TTjl i^-'^i^^ 
However, for a — i: /2 only g^^a is the same in Fi and 
F2. In Fi layer, g„„ is x- and cr-dependent and cannot 
be expressed in a compact form. 

The energy of Andreev states is determined by the 
poles of the Green functions. In contrast to the case 
of collinear magnetizations, the spin splitting of Andreev 
states is absent for orthogonal magnetizations, Eq. ([TU|. 
where the long-range spin-triplet proximity effect is max- 
imally pronounced. 

Real part of the Green functions as a function of en- 
ergy e inside the superconducting gap, — Aq < e < Aq, 
reduces to four delta functions with ^-dependent position 
for collinear magnetizations, and two delta functions for 
orthogonal magentizations. For planar SF1F2S junctions 
these contributions should be numerically summed over 
9. Note that the electronic DOS in ferromagnets is inde- 
pendent of X only in the clean limit. For a finite electronic 



mean free path, DOS is calculated in the middle of F2 
ferromagnetic layer. 



The above results for DOS in the clean limit, obtained 
in the stepwise approximation for the pair potential, are 
not altered qualitatively in numerical self-consistent cal- 
culations. Numerical computation for both the clean 
limit and for moderately diffusive ferromagnets is car- 
ried out using the collocation method: Eq. ^ is solved 
iteratively together with the self-consistency condition, 
Eq. ([S]). Iterations are performed until self-consistency 
is reached, starting from the stepwise approximation for 
the pair potential. For a finite electron mean free path 
in ferromagnets, the iterative procedure starts from the 
clean limit. We choose the appropriate boundary condi- 
tions in superconductors at the distance exceeding 2^^ 
from the SF interfaces. These boundary conditions are 
determined by eliminating the unknown constants from 
the analytical solutions in stepwise approximation. To 
reach self-consistency with good accuracy, starting from 
the stepwise A, five to ten iterative steps were sufficient. 
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FIG. 2. (Color online) Normalized DOS, iV(e)/iVo(0), in ferro- 
magnets for an asymmetric SF1F2S junction with di — lOkp^ 
and d2 = 990kp^ , h\ = h2 = O.IEf- Superconductor is char- 
acterized by Ao(0)/i5F = 10~^; the phase difference (f> — 0. 
The results of non-self-consistent calculations in the clean 
limit, Z — >■ 00, are shown for three values of the relative angle 
between magnetizations: a — (dashed curve), 7r/2 (solid 
curve), and tt (dash-dotted curve). 



III. RESULTS 

We illustrate our results on SF1F2S planar junctions 
with relatively weak ferromagnets, h/Ep = 0.1, and the 
ferromagnetic exchange length S^p = 20fc^^. Supercon- 
ductors are characterized by the bulk pair potential at 
zero temperature Ao(0)/£'f = lO^'^, which corresponds 
to the superconducting coherence length ^s(O) = 636A:^^. 
We assume that all interfaces are fully transparent and 
the Fermi wave numbers in all metals are equal {k^^ ^ 

lA). 

Detailed analysis is given for a highly asymmetric junc- 
tion (di = lOfc^^ and d2 = 990fc^^), Figs. [H and for a 
symmetric junction (di = d2 ~ 500fcj^^), Figs. [S]and[Sl 
In these examples: hi = h2 = O.lEp and T/Tc = 0.1. In 
Figs. [21 [31 and[5l the results of non-self-consistent calcu- 
lations for DOS and Joscphson currents are shown in the 
clean limit {I — >■ 00). The self-consistent numerical cal- 
culations for electronic DOS in the clean limit {I — > 00) 
and for moderate disorder {I = 200k p^) in ferromagnets 
are shown in Figs. [31 and [51 

In highly asymmetric SF1F2S junctions the variation 
of DOS for finite electronic mean free path is less pro- 
nounced but qualitatively similar to the clean limit re- 
sults, see Fig. [H This is not the case for symmetric 
junctions where the long-range spin-triplet proximity ef- 
fect is weak. Fig. [HI 

Influence of the ferromagnet is predominantly deter- 
mined by the parameter O = {h/ Ep)kpd, and the re- 
sults shown here remain the same for a large range of the 
junction parameters where Oi 1 in Fi and O2 100 
in F2 layers. Note that Oi ^ 1 is the optimal choice for 



FIG. 3. (Color online) The Josephson current-phase relation 
I{(j>) for an asymmetric SF1F2S junction with di — Wkp^ 
and d2 = 990k~^ , hi ^ h2 = O.lEp, I 00, T/Tc = 0.1, and 
for three values of the relative angle between magnetizations: 
a — (dashed curve), 7r/2 (solid curve), and n (dash-dotted 
curve). Inset: Magnified curves for a = and a — tt. 
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FIG. 4. Normalized DOS, iV(e)/A''o(0), in ferromagnets for 
an asymmetric SF1F2S junction with d\ — lOfc^^ and d2 = 
^ ^ — ^ _ rv 1 c — _ 7r/2. Superconductor 

^; the phase difference 



990fc^\ hi = h2 ^ O.IEf, for a 



is characterized by Ao(0)/-Bf — 10" 
</> = 0. Results of self-consistent numerical solutions in the 
clean limit, I — >■ 00 (thin solid curve) and for moderately 
diffusive ferromagnets, / — 200k , x = d2/2 (thick solid 
curve) are shown. For comparison, the non-self-consistent 
solution in the clean limit (dotted curve) is also shown. 



long-range triplet Josephson effect.^'' For experimental 
realization of asymmetric junctions it is more suitable 
to use weak and strong ferromagnets {hi <^ /i2) with 
small and comparable thickness (di ~ c?2) to avoid mul- 
tidomain magnetic structures and destructive influence 
of the orbital effect (vortices) 41 

We have verified by additional calculations that the 
long-range triplet proximity effect in highly asymmet- 
ric SF1F2S junctions is rather robust. Small variations 
in the layers thickness and in the exchange energy of 
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FIG. 5. (Color online) Normalized DOS, iV(£)/iVo(0), in fer- 
romagnets for a symmetric SF1F2S junction with di = = 
aOOkp^ , /ii = /i2 = 0.1-Ef. Superconductor is characterized 
by Ao(0)/-Ep = 10"''; the phase difference = 0. Results 
of non-self-consistent calculations in the clean limit, Z — J- 00 
are shown for three values of the relative angle between mag- 
netizations: a = (dashed curve), ■k/2 (solid curve), and tt 
(dash-dotted curve). 
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FIG. 6. Normalized DOS, iV(e)/iVo(0), in ferromagnets for 
a symmetric SF1F2S junction with di — d2 — 500fc]^^, hi — 
/i2 = 0.1-Ef, for a = 7v/2. Superconductor is characterized 
by Ao(0)/£'f = 10~^; the phase difference (j) ~ 0. Results of 
self-consistent numerical solutions in the clean limit, / — > 00 
(thin solid curve) and for moderately diffusive ferromagnets, 
I — 200fc^^, X = d2/2 (thick solid curve) are shown. For 
comparison, the non-self-consistent solution in the clean limit 
(dotted curve) is also shown. 



ferromagnets, as well as moderate disorder, do not af- 
fect distinctive feature of DOS. Likewise, the long-ranged 
second harmonic in /(</)) at low temperatures remains 
dominant 

In all calculations we have assumed transparent SF 
interfaces and clean or moderately disordered ferromag- 



nets. Finite transparency, as well as large disorder, 
strongly suppress the Josephson current and the elec- 
tronic DOS approaches its normal state value. 

In symmetric SFFS junctions with equal ferromagnetic 
layers no substantial impact of long-range spin-triplet su- 
perconducting correlations on the Josephson current has 
been found previously both in the ballistic and diffusive 
regimes ji^"— It can be seen from Figs. [5] and H] that the 
electronic DOS for a = tt/2 is much less pronounced 
than for the corresponding (equal strength and total 
thickness of ferromagnets) asymmetric junction. Figs. [2] 
and [31 Moreover, for moderate disorder in ferromagnts, 
I — 200kp^, the proximity effect is practically absent, 
and DOS equals its normal value (Fig. In symmetric 
SFFS junctions the Josephson critical current monoton- 
ically increases with angle between the magnetizations, 
except in the vicinity of — tt transitions.-- In the clean 
limit for a = tt, the influence of opposite magnetiza- 
tions in F layers practically cancels out and the proxim- 
ity effect is much stronger, see Fig. [S] Both DOS and 
the Josephson current-phase relation are the same as for 
the corresponding SNS junction, where N stands for a 
normal nonmagnetic metal (hi = /12 = 0)- Recently, 
substantially enhanced supercurrent has been observed 
in Josephson junctions with ferromagnetic Fe/Cr/Fe tri- 
layer in the antiparallel configuration, without long-range 
spin-triplet correlations!^ 



IV. CONCLUSION 

We have studied the spin-triplet proximity effect in 
SF1F2S planar junctions made of conventional supercon- 
ductors and two monodomain ferromagnetic layers with 
arbitrary thickness, strength, and angle between in-plane 
magnetizations. We have derived analytical expressions 
for the quasiclassical normal Green functions in the clean 
limit and the stepwise pair potential approximation (non 
self-consistent). We have calculated the electronic den- 
sity of states in ferromagnets, and the Josephson current. 
In addition, results of the numerical self-consistent calcu- 
lations for the clean and moderately diffusive ferromag- 
nets are shown for comparison. 

The second harmonic in the Josephson current-phase 
relation is dominant for highly asymmetric SF1F2S junc- 
tions composed of particularly thin (weak) and thick 
(strong) ferromagnetic layers with noncollinear magne- 
tizations. This is a manifestation of the long-range spin- 
triplet proximity effect at low temperatures, related to 
the phase coherent transport of two Cooper pairs. We 
find that the proximity effect is also accompanied by 
distinctive variation of electronic DOS in ferromagnetic 
layers as a function of energy near the Fermi surface. 
This variation is less pronounced in moderately diffusive 
case, but still detectable by tunneling spectroscopy as a 
signature of the long-range proximity effect. The self- 
consistency and finite (but large) electronic mean free 
path do not alter the analytical results qualitatively. In 
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contrast, in symmetric junctions with moderately diffu- 
sive and thick ferromagnets the long-range proximity ef- 
fect is absent and DOS is equal to its normal metal value. 

In summary, fully developed long-range proximity ef- 
fect occurs in highly asymmetric SF1F2S junctions with 
orthogonal magentizations. Dominant second harmonic 
in the Josephson current-phase relation, as well as a dis- 
tinctive variation of DOS in ferromagnetic layers with 
the angle between magnetizations, should be experimen- 
tally observable for relatively small interface roughness 
and relatively clean ferromagnetic layers at low temper- 
atures. 
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VI. APPENDIX 

In the clean limit, ^ cxo, solutions of Eq. ([T|) in the 
step- wise approximation, Eq. ([8]), for the left supercon- 
ductor [x < — c?i) can be written in the usual form for 
normal Green functions 



and for anomalous Green functions 



/ti(^) 



2uJn + hVxKs 

iAo 2iAo 



'2fc 



Qn 2uJn + hVxKs 



rin 2uJn + hVxKs 



with 



2W„ + hv^Ks 



HvfI cos(^)| ' 



(12) 

(13) 
(14) 

(15) 



(16) 
(17) 
(18) 
(19) 

(20) 



and rtn = \/oJn + |Ap. For the right superconductor 
(x > £^2), the solutions retain the same form with Ks — > 
—Kg, 4> — ^ with a new set of constants D'^^ ...,D'^. 

Solutions for the Green functions in the left ferromag- 
netic layer Fi, —di < x < 0, for three values of angle 
ai — 0, 7r/2, TT (we have chosen a2 — 0) can be written 
as follows. For collinear magnetizations ai = and tt, 



and 



gti(x) = if2e±'^'="% 
g^^(x)=if3e™, 

fnix) = Cle^''^ 
fnix) = C2e''=±% 
Mx) = Cse'-^^, 



(21) 
(22) 
(23) 
(24) 

(25) 
(26) 
(27) 
(28) 



where the upper (lower) sign corresponds to ai = (ai 
7r). For orthogonal magnetizations, ai = 7r/2, 



and 



9ni^) 

9lt{x) 
9ui^) 

fuix) -■ 



K4 + 1X26""""= 



Here, 



kq = 2hi/hvx, 

K = 2iuin/hVx, 



(29) 
(30) 

(31) 
(32) 

(33) 
(34) 
(35) 
(36) 

(37) 
(38) 
(39) 



For the right ferromagnetic layer F2, < a; < (i2, 0^2 = 
0, the solutions are the same as for the layer Fi in the 
case ai — 0, with hi — > h2 and a new set of constants 
K[, . . . , K'^, C'l, . . . , C4. 

The complete solution has 24 unknown coefficients: 
4 -|- 4 in superconducting electrodes, and 8 + 8 in fer- 
romagnetic layers. Continuity of normal and anomalous 
Green functions at the three interfaces provide the nec- 
essary set of 24 equations. 



1 F. S. Bergeret, A. F. Volkov, and K. B. Efetov, Phys. Rev. 

Lett. 86, 4096 (2001). 
^ A. F. Volkov, F. S. Bergeret, and K. B. Efetov, Phys. Rev. 



Lett. 90, 117006 (2003). 
^ F. S. Bergeret, A. F. Volkov, and K. B. Efetov, Rev. Mod. 
Phys. 77, 1321 (2005). 



7 



* 1. Sosnin, H. Cho, V. T. Petrashov, and A. F. Volkov, 
Phys. Rev. Lett. 96, 157002 (2006). 

5 R. S. Keizer, S. T. B. Goennenwein, T. M. Klapwijk, 

G. Miao, G. Xiao, and A. Gupta, Nature 439, 825 (2006). 
® M. S. Anwar, F. Czeschka, M. Hesselberth, M. Porcu, and 

J. Aarts, Phys. Rev. B 82, 100501 (2010). 
^ J. Wang, M. Singh, M. Tian, N. Kumar, B. Liu, C. Shi, 
J. K. Jain, N. Samarth, T. E. MaUouk, and M. H. W. 
Chan, Nature Phys. 6, 389 (2010). 

* D. Sprungmann, K. Westerholt, H. Zabel, M. Weides, and 

H. Kohlstedt, Phys. Rev. B 82, 060505(R) (2010). 

^ T. S. Khaire, M. A. Khasawneh, W. P. Pratt, and N. O. 

Birge, Phys. Rev. Lett. 104, 137002 (2010). 
^° J. W. A. Robinson, J. D. S. Witt, and M. G. Blamire, 

Science 329, 59 (2010). 
" V. Braude and Y. V. Nazarov, Phys. Rev. Lett. 98, 077003 

(2007). 

M. Houzet and A. I. Buzdin, Phys. Rev. B 76, 060504 

(2007) . 

A. F. Volkov and K. B. Efetov, Phys. Rev. B 81, 144522 
(2010). 

Z. Pajovic, M. Bozovic, Z. Radovic, J. Cayssol, and 

A. Buzdin, Phys. Rev. B 74, 184509 (2006). 

B. Crouzy, S. ToUis, and D. A. Ivanov, Phys. Rev. B 75, 
054503 (2007). 

^® L. Trifunovic, Z. Popovic, and Z. Radovic, Phys. Rev. B 
84, 064511 (2011). 

L. Trifunovic, Phys. Rev. Lett. 107, 047001 (2011). 
Z. Radovic, L. Dobrosavljevic-Grujic, and B. Vujicic, 
Phys. Rev. B 63, 214512 (2001). 
A. L Buzdin, Rev. Mod. Phys. 77, 935 (2005). 
^° M. Eschrig, J. Kopu, J. C. Cuevas, and G. Schon, Phys. 
Rev. Lett. 90, 137003 (2003). 

M. Eschrig and T. Lofwander, Nature Phys. 4, 138 (2008). 
^2 A. F. Volkov and K. B. Efetov, Phys. Rev. B 78, 024519 

(2008) . 

L. Trifunovic and Z. Radovic, Phys. Rev. B 82, 020505 

(2010) . 

2* M. Alidoust and J. Linder, Phys. Rev. B 82, 224504 (2010). 
M. Duckheim and P. W. Brouwer, Phys. Rev. B 83, 054513 

(2011) . 

A. Cottet, Phys. Rev. Lett. 107, 177001 (2011). 

T. Yokoyama, Y. Tanaka, and A. A. Golubov, Phys. Rev. 



B 75, 134510 (2007). 

T. Yokoyama and Y. Tserkovnyak, Phys. Rev. B 80, 
104416 (2009). 

S. Pilgram, W. Belzig, and C. Bruder, Phys. Rev. B 62, 
12462 (2000). 

^° F. Perez- Willard, J. C. Cuevas, C. Siirgers, P. Pfundstein, 
J. Kopu, M. Eschrig, and H. v. Lohneysen, Phys. Rev. B 
69, 140502 (2004). 

Y. Asano, Y. Tanaka, and A. A. Golubov, Phys. Rev. 
Lett. 98, 107002 (2007). 

K. Halterman, P. H. Barsic, and O. T. Vails, Phys. Rev. 
Lett. 99, 127002 (2007). 

K. Halterman, O. T. Vails, and P. H. Barsic, Phys. Rev. 
B 77, 174511 (2008). 

Y. Fominov, A. Golubov, and M. Kupriyanov, JETP Lett. 

77, 510 (2003). 

C. You, Y. B. Bazaliy, J. Y. Gu, S. Oh, L. M. Litvak, and 
S. D. Bader, Phys. Rev. B 70, 014505 (2004). 
T. Lofwander, T. Champel, J. Durst, and M. Eschrig, 
Phys. Rev. Lett. 95, 187003 (2005). 

Y. S. Barash, I. V. Bobkova, and T. Kopp, Phys. Rev. B 
66, 140503 (2002). 

Y. V. Fominov, A. F. Volkov, and K. B. Efetov, Phys. 
Rev. B 75, 104509 (2007). 

L B. Sperstad, J. Linder, and A. Sudb0, Phys. Rev. B 78, 
104509 (2008). 

40 rp Y Karminskaya, A. A. Golubov, M. Y. Kupriyanov, 
and A. S. Sidorenko, Phys. Rev. B 79, 214509 (2009). 

'^^ Y. M. Blanter and F. W. J. Hekking, Phys. Rev. B 69, 
024525 (2004). 

J. Linder, M. Zareyan, and A. Sudb0, Phys. Rev. B 79, 
064514 (2009). 

G. Eilenberger, Z. Phys. 190, 142 (1966). 
Our matrix g and Eq. ^ are equivalent up to an unitary 
transformation to the definitions used in Ref. [^. 
B. Miihlschlegel, Z. Phys. 155, 313 (1959). 
'^^ F. Konschelle, J. Cayssol, and A. I. Buzdin, Phys. Rev. B 

78, 134505 (2008). 

L F. Lyuksyutov and V. L. Pokrovsky, Adv. Phys. 54, 67 
(2005). 

J. W. A. Robinson, G. B. Halasz, A. I. Buzdin, and M. G. 
Blamire, Phys. Rev. Lett. 104, 207001 (2010). 



